Power-law tails from multiplicative noise 
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Abstract We show that the well-known linear Langevin equation, modeling the Brownian mo- 
tion and leading to a Gaussian stationary distribution of the corresponding Fokker-Planck equation, 
is changed by the smallest multiplicative noise. This leads to a power-law tail of the distribution for 
large enough momenta. At finite ratio of the correlation strength for the multiplicative and additive 
noise the stationary energy distribution becomes exactly the Tsallis distribution. 



Power-law tails are present in numerous distributions 
studied in physics or elsewhere when dealing with com- 
plex systems They are of generic interest, regarded as 
to signal long range order, non- vanishing correlations or 
scale invariance in complex systems with strong dynam- 
ics, which's details are mostly unknown. These are often 
contrasted to the traditional statistical system, showing 
the Gibbs distribution (exp(— E/T)) in energy, which 
is Gaussian in the momenta of free, massive particles 
(exp(— p 2 /2mT)). The latter is considered as the generic 
case for thermal equilibrium of non-correlated or short- 
range correlated systems. This concept has been carried 
far beyond of its original field describing monoatomic 
ideal gas (Maxwcll-Boltzmann statistics), by applying 
the Gibbs distribution in thermal equilibrium to areas 
such as particle physics and field theory. It serves as 
starting point of high-temperature field theory calcula- 
tions both with analytical and numerical methods. Lat- 
tice gauge theory is based on the formal similarity be- 
tween Euclidian path integrals and the canonical parti- 
tion sum. 

A very simple and elegant, microdynamical explanation 
for the Maxwcll-Boltzmann statistics is offered by the 
Langevin equation, describing a free particle moving un- 
der the influence of a deterministic damping force and 
a stochastic drive 0- The latter accelerates the parti- 
cle in a short time, changing its momenta randomly and 
uncorrelated. The stationary solution of this stochastic 
equation follows the Gaussian statistics, compatible to 
Gibbs' principle. It seems that in many statistical con- 
siderations of complex physical models from that on it is 
tacitly assumed that the presence of this additive noise 
is a dominant effect: the equilibrium distribution follows 
Gibbs' formula (the Gaussian distribution in momentum 
for a free, massive particle). Since the harmonic oscillator 
is just the extension of this free motion Langevin equa- 
tion into the phase space, also for free quantum systems 
the above picture is generally accepted. 

We shared this expectations, thinking that any non- 
Gaussian (or in the energy non-exponential) distribution, 
especially the power-law tail observed in many phenom- 
ena including particle spectra in high energy physics, may 



only come from non-thermalizcd or in an other way non- 
equilibrium situation. In this note we would like to share 
our deep astonishment about that this is not so: we found 
out that treating the damping constant in the Langevin 
equation also stochastically (considering this way both 
multiplicative and additive noise) the stationary distri- 
bution is in general non-Gaussian. Moreover, above a 
certain momentum, depending on the strength (i.e. self- 
correlation width) of the multiplicative noise, the station- 
ary solution goes over into a power-law. 

This finding seems to have manifold consequences. Even 
in equilibrium, even averaging over a huge number of 
elementary events, measurements on quantum systems, 
such as particles or fields, are bound to find a power-law 
tail irrespective if the underlying dynamical system had 
equilibrated itself long enough in units of a characteristic 
time. Power-law tails in pion spectra, found experimen- 
tally in high energy e + e~, pp or heavy- ion collisions Q, in 
the view of this simple mathematical result may reflect 
an already stationary distribution. Power-law tails found 
in other areas, in particular as properties of average cor- 
relations (cf. time series in stock market), may also be a 
non-temporary, long term effect. 

In this paper we derive a generic stationary distribution 
for the Langevin-type equation with both additive and 
multiplicative noise. Conform to the original assump- 
tions both noise terms are white (Dirac delta correlated 
in time), but they may show cross-correlations. This situ- 
ation may stem naturally from the widespread treatment 
of field theoretical operator equations, where the fields 
are split to a large, "classical" expectation value and to 
a noisy quantum or thermal fluctuation part: the differ- 
ent noise terms do have a common origin, so it is natural 
to consider cross-correlations among them £| . Often the 
overdamped approximation is applied for studies of plas- 
mas, solving then effectively a differential equation first 
order in time, instead of second orderj^ . We think there- 
fore that the analytically solvable case, we present in this 
paper, may carry a quite general lesson. 

We present the solution of the Langevin equation with 
both additive and multiplicative noise terms applying the 
classical method of Wang and Uhlenbeck0|. Our starting 



2 



point is the linear equation: 

P + IP = & 



(1) 



where now both £ and 7 are stochastic variables. They 
both may have a non-zero mean value (motivated by pos- 
sible field theory applications), 



W))=F, (j(t))=G, 



(2) 



and show white-noise (i.e. extremely short term) corre- 
lations: 

(7(t)7(t'))-(7(t))( 7 (t')) = 2C6(t-t'), 

mw))-m)(7(t'))=2BS(t-t'), (3) 

This problem can be solved analytically. We will deter- 
mine the time dependence of the distribution of p values, 
denoted by f(p, t). In this notation f(po, t)dp is the prob- 
ability that after time t the variable p has the value in 
the range [po,Po + dp]. We rewrite Q as a difference 
equation 



t+dt 



p(t + dt) = p{t) 



dt' (e(0-7(0K0), (4) 



If p(t) is a smooth function of time we can replace pit 1 ) 
either by p(t) or p(t + dt), or any value in between. We 
choose here the Ito prescription [jj, which uses p(t) under 
the integral in the — > limit. In order to simplify nota- 
tion we write the integral term as dt (x) , with x denoting 
the general integrand. Now we can write down a Fokker- 
Planck equation for the distribution: the probability to 
have the value pit + dt) at t + dt is the probability that 
we have the value pit) at time t and noise values £ and 
7 that just satisfy 10} : 

f(p,t + dt)= fd^d 7 P^,7)fip-dt{0+pdt{ 7 )). (5) 

Unfortunately this form is not appropriate to directly 
create differential equation as dt — > 0. Instead we follow 
the method of Wang and Uhlenbeck: we have a trial 
function Rip) that is smooth enough and we compute 
the expectation value of R (averaged over the noise) as a 
function of time 



(R(t)) = jdpR(p)fip,t). (6) 
Applying this form to J5| we have 

JdpR(p)f(p,t + dt) = 
Jdtdr/P(t,i) JdpRip + dt (0-pdt <7»/(p,*).(7) 



By Taylor expanding Rip) and integrating over the noise 
distribution we get 

(Rip + dt (O-Pdt (7))) = 

(Rip) + dtR\p) K x ip) + dtR"ip) K 2 (p) + Q(dt 2 ))(8) 
with 



A'i = F-Gp, K 2 = D- 2Bp + Cp 2 



(9) 



in the present case. This leads to the following general 
Fokkcr-Planck equation 



df 
dt 



diK\f) + d\K 2 f) 



dp dp 2 
The stationary solution satisfies 
d 



dp 



(A 2 /) = Ai /, 



which is analytically solvable. It leads to 



fip) = f(0)^eMHp)) 



with 



Lip) = ( P dq 
Jo 



KM 
A 2 (g)' 



(10) 



(11) 



(12) 



(13) 



In the case of two noises correlated the way given in eq. © 
we arrive at the following logarithm of the stationary 
distribution: 



In 



/(0) 
with 



= - 1 
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2C 
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K 2 jp) 
D 



— — atn 



dp 



D - Bp 



d = \JCD - B 2 



and 



= G 



B 
C 



F. 



(14) 



(15) 



Here 'atn' denotes the inverse tangent function, not al- 
ways taking the first principle value, but rather contin- 
uing at p > D/ B smoothly. Considering physical appli- 
cations the noise correlation values, D, C and B build 
a positive semidefinite matrix. This ensures that C and 
D are non-negative values, and the determinant $ is real 
and also non-negative. The same applies for the function 
K 2 ij>) occurring under the logarithm. In this context zero 
values are limiting cases and the stability of the station- 
ary solution (|14|) against choosing a small positive value 
has to be investigated. 

The lesson for physical applications lies in the analysis of 
different limiting cases. First we consider the traditional 
case: C = B = 0, not allowing for any noise (fluctuation) 
in the multiplicative factor 7 (damping constant). This 
leads back to the familiar Gauss distribution: 



M = /(0)e 



- — 7)' 

2D f 



-I' 



(16) 
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FIG. 1: Comparison of the generic stationary distributions 
of the Langevin-type equation with both additive and multi- 
plicative noise. Limiting cases, as the Gauss-, the Gamma- 
and the Power-Law distribution are labeled correspondingly. 
The label 'Zero' refers to the case with •& = 0. 



with an eventual shift in the mean momentum for a non- 
zero mean driving force (£) = F. Another limiting case is 
that of the purely multiplicative noise with D = 0, B = 
(G.Wilk|| has applied it for the heat conduction equa- 
tion and obtained a Gamma distribution for the inverse 
temperature 1/T). Now the stationary solution becomes 



/(p)=/(0)fT 



G/C e c | p | 



(17) 



a Gamma distribution in 1/p. For large p this ap- 
proaches a pure power-law. It is particularly interest- 
ing to investigate the mathematically degenerate case of 
t) = 0. Now K2 (p) reaches zero at the critical momen- 
tum, p c = \J D/C, leading to zero probability in the 
stationary distribution f(y/D/C) = 0. This occurs as a 
"limiting momentum" in the physical distribution. 
The above limiting cases rely on an expansion of the 
generic solution. They are valid only in a limited range of 
momenta: the Gaussian solution for small, the Gamma 
distribution for large argument of the inverse tangent 
function. Correspondingly the widely beloved Gaussian 
distribution can be a good approximation to the sta- 
tionary solution only for momenta p <C y/D/C. (For 
C = strictly, of course this is 'the' solution for any fi- 
nite momentum.) Generally the small argument of the 
inverse tangent is fulfilled for p <C ^xs" This result also 
means that for the smallest fluctuation in the multiplica- 
tive factor the stationary Gauss distribution develops a 
power-law tail. The power in this tail, p~ 2v , is given by 
v = 1 + . In order to offer a visual insight into the na- 
ture of the generic stationary solution we show stationary 
spectra for different parameters (cf.FigJTJ). 
In the case B = (no cross correlation between the 
noises), and F = (no drift term due to the additive 
noise), the exact stationary distribution is the Tsallis 
distribution 9] . In order to achieve this result one uses 
the energy of the free particle, E = p 2 /2m as the distri- 



bution variable and ea. (|14f) . We get 



l + (?-l)f V 



(18) 



where the parameters of the Tsallis distribution are given 
by 



T 



D 

mG ' 



1 



2C 
G ■ 



(19) 



Again for C = (only additive noise) q = 1 and the 
Tsallis distribution goes over into the Gibbs distribution, 



f(E) = foeM-E/T). 



(20) 



Tsallis and others have worked out a thermodynamical 
framework offering the distribution (|18|) as the canonical 
distribution. This approach, the non-extensive thermo- 
dynamics, however, is based on a non-extensive entropy 
measure. This eventually unwanted property is not fatal: 
the distribution (|18|) can also be obtained based on the 
extensive Renyi entropy [loj. We note by passing this 
point that the presence of the two uncorrelated noise 
and the corresponding Fokker-Planck equation also can 
be obtained from an inhomogeneous diffusion coefficient. 
Instead of JD one may consider 



p + Gp=(D + Cp 



2N1/2 



(21) 



With a single noise n, normalized to unity: 

<»/(*)> = 0, {r,{t)n{t>))=28{t-U). (22) 

This is a particular case of a more general field-dependent 
noise considered by Arnold, Son and Yaffe in the context 
of non-abelian plasmas|llj. 

The mechanism outlined in the present article may work 
for non-abelian gauge theories, as well. The basic glu- 
onic field, described by a vector potential A, satisfies the 
equation 



A + a A = D x (D x A), 



(23) 



with D being the gauge-covariant derivative and a the 
color conductivity factor. We consider an overdamped 
dynamics, when the second time derivative is ignored. 
Driving this to the extreme we single out in a Fourier 
expansion the zero mode, Aq and consider the following 
effective equation 

o-A = J2(k-gA k )(q-gA q )A_ k _ q . (24) 

k,q 

Treating the the k = and q = contributions on the 
right hand side separately we arrive at 



^0 - — Al = -jA Q + 
a 



(25) 
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FIG. 2: 

with 

7 = -2~^2(k - gA k )A- k , 
fe>o 

£ = - E (26) 

fc>0,<?>0 

Ignoring the classical ~ A\ contribution for the follow- 
ing discussion, one realizes that 7 and £ do contain hard 
Fourier component contributions. In the Langevin equa- 
tion with both additive and multiplicative noise these are 
regarded as noisy, fast-fluctuating quantities. 
The full quantum field theoretical treatment of these fac- 
tors is rather involved. In order to gain estimates over 
the parameters of a possible stationary distribution of 
Aq values, averages and correlations of 7 and £ have to 
be obtained. This process can be facilitated by using a 
graphical notation: external legs stand for the zero mode 
and internal lines for the hard modes, (cf. Fig[2|>. For 
no spontaneous symmetry breaking, F = and B = 
immediately follows, so the stationary distribution is a 
Tsallis-distribution for Aq and eventually in the soft en- 
ergy E ~ VK 2 A 2 = A 2 /(2m), where A is the underly- 
ing energy scale. A qualitative, order of magnitude es- 
timate for the parameters of the Tsallis distribution can 
be given as: T — g 2 A 4 /(mcr 2 ) and q — 1 ~ g 2 A 2 /a 2 . 
The starting point of the power law in ca. (|18fl is at 
E c = T/(q - 1) ~ A 2 /m. 

Finally we would like to check whether quantitative es- 
timates give reliable results. In high energy particle 
physics experiments the transverse momentum distri- 
bution has been investigated for long. From Gaussian 
fits to the parton distribution one conjectures a ra- 
tio D/G = (pi) w 1 - 1.5 GeV 2 . On the other 
hand power-law tails at high transverse momenta make 
a value of v « 5.8±0.5 realistic. This fixes the ratio to 
G/C « 9.6±1 and the Tsallis index to q « 1.2±0.03. The 



critical transverse momentum, beyond which the power- 
law dominates the familiar Gauss distribution, can be 
calculated from this to be p c « 3 — 4 GeV. Compared to 
experimental spectraQ this is a quite realistic estimate. 

In conclusion we have shown that the smallest multi- 
plicative white-noise related to the classically determin- 
istic damping constant in the Langevin equation leads 
to a stationary distribution of particle momenta, which 
ends in a power-law tail at high values. This result seems 
to undermine the well-established thermal approaches to 
phenomcnological and field theoretical studies in parti- 
cle physics, where the presence of a multiplicative noise 
is not less probable than the presence of an additive one 
in any simplification (linearization) of the underlying mi- 
crodynamical problem. This approach on the other hand 
offers a new way to deal with the interpretation of power- 
law tails occurring in experimental findings, as well as it 
animates to seek new methods in thermal field theory 
exceeding the traditional thermodynamical approach. 
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